perm filename ELLIP2.F4[SAB,LCS] blob sn#349440 filedate 1978-04-16 generic text, type C, neo UTF8
COMMENT āŠ—   VALID 00002 PAGES
C REC  PAGE   DESCRIPTION
C00001 00001
C00002 00002		SUBROUTINE ELLIP2(A,B,XC,YC,PSI)  
C00003 ENDMK
CāŠ—;
	SUBROUTINE ELLIP2(A,B,XC,YC,PSI)  
	N=IFIX(A*50.)+12
	IF(B.GT.A)N=IFIX(B*50.)+12
	PHI= (6.2831853/360.)*PSI
	X=A
	Y=0.0
	THETA=6.2831853/FLOAT(N)
	C1=COS(THETA)
	C4=A/B
	C6=SIN(THETA)
	C2=C6*C4
	C3=C6/C4
	C=COS(PHI)
	S=SIN(PHI)
	CALL PLOT(A*C+XC,A*S+YC,3)
	DO 10 I=1,N
 	U=X*C1-Y*C2
	V=X*C3+Y*C1
	X=U
	Y=V
	X1=X*C-Y*S
	Y1=X*S+Y*C
	CALL PLOT(X1+XC,Y1+YC,2)
 10	CONTINUE
	RETURN
	END